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Chemical control of the spontaneous motion of a reactive oil droplet moving on a glass substrate 
under an aqueous phase is reported. Experimental results show that the self-motion of an oil droplet 
is confined on an acid-treated glass surface. The transient behavior of oil-droplet motion is also 
observed with a high-speed video camera. A mathematical model that incorporates the effect of the 
glass surface charge is built based on the experimental observation of oil-droplet motion. A numerical 
simulation of this mathematical model reproduced the essential features concerning confinement of 
oil droplet motion within a certain chemical territory, and also its transient behavior. Our results 
may shed light on physical aspects of reactive spreading and a chemotaxis in living things. 
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I. INTRODUCTION 



The spontaneous motion of reactive droplets has attracted considerable attention in relation to energy transduc- 
tion by living organisms, i.e. chemomechanical energy transduction. In 1978, the motion of a droplet driven by a 
surface-tension gradient was predicted by Greenspan to explain some aspects of cell cleaving from a physicochemical 

n 

viewpoint 1J. Various kinds of droplet motion driven by a surface-tension gradient have been reported since, where the 



ay 

or a droplet driven by a sel: 
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asymmetry in a surface condition. Spontaneous 



surface-tension gradient was given by an initial |2j,|3| or an external 

motion in reactive spreading, or a droplet driven by a self-produced surface-tension gradient, has also been noted in 

for example, a droplet of surfactant on a gold surface 



various experimental systems 

Q 

alloy melting |6j, and a water droplet containing silane on a glass substrate 



^y 



To interpret these phenomena, 

extensive theoretical studies on reactive droplets have Ikvii performed with I he lubrication approximation |lJ. lfil lfij . 
An oil-water system composed of an organic phase with potassium iodide and iodine and an aqueous phase containing 



stearyl trimethyl ammonium chloride (STAC) exhibits self-agitation at the oil- water interface, accompanied by spatio- 



temporal instability of interfacial tension [17 




droplet in an oil-water system exhibits chemo-sensitivity 



21 

23 



2^. It has also been found that the motion of am oil 
Based on the measurement of the electrical potential 



at the oil- water interface in a similar system, the nature of the electrical fluctuation/oscillation has been shown to 
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strongly depend on the chemical properties and the concentration |24| . However, the detailed mechanism of this effect 
of chemical substances has not yet been clarified. Recently, it was proposed that an oil-water system also shows 
reactive spreading on a glass surface with recovery of the surface condition 

In the present study, we experimentally observed the chemo-sensitive motion of an oil droplet on a glass substrate 
that had been partially pretreated with acid. We found that an oil droplet exhibits various behaviors, such as turning- 
back motion, stopping at a certain position, or slowing down around regions treated with acid. In addition, we also 
noted that the velocity and shape of the oil droplet show damped oscillation while the droplet moves continuously. 
However, detailed observations show that the motion of an oil droplet strongly depends on slight changes in the 
surface condition of the substrate. In terms of the argument given in our previous study 2fjj we expect that the 
chemo-sensitivity of oil-droplet motion is caused by the condition of the glass surface affected by acid, which inhibits 
the aggregation of stearyl trimethyl ammonium ion (STA + ion) on the glass surface. To confirm this hypothesis, we 
idealize an oil droplet as self-propelling spring-beads whose driving force is given by the surface-tension gradient. This 
model successfully reproduces the observed experimental results. Furthermore, a numerical simulation of the model 
reveals that the selection among the turning-back motion, stopping at a certain position, and slowing down can be 
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FIG. 1: Schematic representation of the experimental setup (a) on a Y-shaped path on a glass surface, and (b) on an acid- 
treated narrow glass substrate. Prior to the experiment, both ends of the narrow glass substrate (ft; 15 mm) were treated with 
acid. 



explained in terms of the effect of the acid on the glass substrate, where minute changes in the acid effect drastically 
change the behavior of the droplet. As has been reported previously, the oil- water system exhibits temporal fluctuation 
in electrochemical potential, which is sensitive to chemical stimuli |24j. On the other hand, the present study shows 
that it may be possible to acquire spatio-temporal information on minute changes in the chemical concentration on 
a surface by noting a drastic change in the characteristic motion of a self-running oil droplet. This result may be 
associated with the mechanism of chemotaxis in biological systems, where living organisms show extreme sensitivity 
to a small spatial gradient in chemicals [27 , . 



II. EXPERIMENTS 

The aqueous phase contained 1 mM STAC and the organic phase was 5 mM iodine solution of nitrobenzene saturated 
with potassium iodide. STAC was prepared by recrystallization using acetone. For glass substrates, we used micro 
slide glass (Matsunami, Osaka; S9111). The glass surface was treated with acid as follows. First, 1 M sulfuric acid 
solution was dabbed on the glass substrate with a cotton swab, and then it was rinsed from the surface with running 
distilled water. Extreme care was taken in rinsing to avoid invading the bare glass surface. 

In the measurement on a Y-shaped path, 100 /il of the oil was placed in an aqueous phase on the pretreated glass 



4 

(a) 




Acid Acid 0.1 sec per frame 




Acid Acid 0.1 sec per frame 



FIG. 2: Experiments on spontaneous droplet motions, (a) Confined motion of an oil droplet. The path of an oil droplet (100 
/il) was enclosed within the acid-treated surface. The left image shows the actual appearance of a substrate. The outside of the 
Y-shaped region was treated with acid. In the middle image, the path of the CM of an oil droplet is shown. The superposition 
of both images is shown on the right. (b),(c) Motion of a 15-/il oil droplet on a narrow acid-treated glass substrate. An image 
from the experiment and a spatio-temporal image of oil-droplet motion are shown. An oil droplet showed (b) shuttling motion 
and (c) intermittent shuttling motion. When an oil droplet did not intrude into the acid-treated region it showed turning-back 
motion, which resulted in shuttling motion. 

substrate [Fig. 1(a)]. A digital video camera (Panasonic; NV-GS100K-K) was used to record the motion of an oil 
droplet at 30 frames per second. 

In the experiments on an acid-treated narrow glass substrate, the glass substrate (1 mm x 52 mm x 3.1 - 3.2 mm) 
was cut from slide glass. Both ends of this substrate were pretreated with acid for 15 mm [Fig.QJb)]. An oil droplet 
of 15-60 fil was placed on the substrate in the same manner as in the Y-shaped path experiment. The motion of an 
oil droplet was recorded by a high-speed video camera (RedLake MASD Inc., San Diego, CA; Motion Scope PCI) at 
125 frames per second. All measurements were carried out at room temperature (20°C ±3°C). The recording of a 
droplet motion was digitized by an image-processing system (Library, Tokyo) . 
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FIG. 3: Experiments on spontaneous droplet motion. Motion of a 15-^1 oil droplet on a narrow acid-treated glass substrate. An 
image from the experiment and a spatio-temporal image of oil-droplet motion are shown. An oil droplet showed (a) stopping 
and (b) slowing. When an oil droplet intruded into the acid-treated region, it showed stopping (a) or slowing (b). 

III. RESULTS 

The spontaneous motion of an oil droplet was found to be confined within a Y-shaped path and a narrow substrate 
[Figs. [21 El- On the Y-shaped path [Fig. Ufa)], motion is clearly limited to within the region that was not treated 

rn 

with acids. On the narrow substrate, oil droplets exhibit quasi 1-dimensional motion 26], i.e. an oil droplet of 15 ^1 
showed turning-back motion that resulted in a shuttling motion [Fig. |2[b),(c)] when the oil droplet could not intrude 
into the acid-treated region. The time trace of oil-droplet velocity v and the apparent length of an oil droplet r in the 
motion shown in Fig.|2b) are plotted in Fig. ^ v and r showed damped oscillation after each returning motion until 
both reached steady values. On the other hand, when an oil droplet of 15 fjl intruded into the acid-treated region, 
the oil droplet stopped and stood still [Fig. Ola)] or slowed down [Fig. O^b)] . The motion of an oil droplet was found 
to be sensitively dependent on the manner of acid treatment. 



IV. DISCUSSION 



It has been suggested by our previous study 26] that the motion of an oil droplet can be explained in terms of 
reactive spreading 15]. STA + ions in the aqueous phase aggregate on the glass surface that has negative charges at 
neutral pH 2^. The organic phase dissolves STA+ ions that have adhered to the glass substrate, and this dissolution 
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FIG. 4: Time trace of the apparent length, r, and velocity, v, of an oil droplet for the motion shown in Fig. |5Jc) 




FIG. 5: Schematic diagram of oil-droplet motion and glass surface conditions. STA + is represented as a hydrophobic bar with 
a hydrophilic head. 



is promoted since STA + and 1% in the organic phase make an ion-pair. The ion-pair dissolves well in non-polar solvent, 
i.e. the organic phase. Therefore, a difference in interfacial energy arises from the difference in the concentration of 
STA + ions on the glass surface between the front and back of an oil droplet [Fig. [S], and the oil droplet is driven by 
Marangoni effect. Thus, the condition of the glass surface on which STA + ions aggregate plays an important role in 
the motion of an oil droplet. 

Nakata and coworkers examined the chcmo-sensitivity of the motion of an oil droplet escaping from acids 
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FIG. 6: Schematic diagram of the model of oil-droplet motion on a glass surface. The lower arrow represents the spatial 
coordinate, where to £o and £ — £o to £ are treated with acid. «o and Wi are the saturated concentration of the STA + ions on 
the glass substrate on bare and acid-treated glass substrate, respectively. The oil droplet is modeled as two lines (perpendicular 
to the x-axis), with a mass density of a, connected to a spring with an elastic constant of /3(ro) (see also Appendix A) 



their experiment, hydrochloric acid was injected in the aqueous phase while an oil droplet was moving. From the 
above discussion, modification of a glass surface by acids significantl y aff ects oil-droplet motion. It has been reported 



3(J,|31|. Thus, the glass surface charge may have 



that the point of zero charge (pHo) of the SiC>2 surface is around 2.2 
been partially reversed to be positive due to the effect of acid since the pH of 1M sulfuric acid is around 0. We can 
consider that this results in fewer adsorption sites for STA + ions on the glass surface. This may explain why an oil 
droplet showed confined motion [Fig. Eta)] and shuttling motion [Fig. Elb)] on an acid-treated glass substrate when 
it was encircled by an acid-treated glass surface. Furthermore, the dependence of droplet motion on the condition of 
the acid-treated region can be attributed to the difference in the number of adsorption sites on the glass surface. To 
explain the various motions of an oil droplet, it seems appropriate to control the surface condition of a glass substrate 
much more accurately. On the other hand, in the experimental process, these different types of the characteristic 
motions were observed even under nearly identical conditions, including the concentration and duration of treatment 
with sulfuric acid. To understand this strong dependence on the condition of the glass surface, we constructed a 
mathematical model, and conducted a numerical simulation of oil-droplet motion on a narrow substrate. 

There may be various ways to theoretically analyze for the spontaneous motion of an oil droplet, such as interface 
dynamics based on Young's equation 1J, Il5| or a hydrodynamic approach which adopts the lubrication approxima- 
tion ^j. Instead of making a detailed model, we propose a simple "beads-spring" model for an oil droplet to abstract 
the essential features of droplet motion. 
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V. MODELING 

Since we want to know the basic mechanism of the motion of an oil droplet on a glass substrate, we make a simple 
one-dimensional model. The x axis is set along the surface of the glass substrate. The droplet proceeds on a line with 
stretching and shrinking. The positions at the both ends of the droplet are set as X\(t) and x 2 (t), where x\ < x 2 . The 
concentration of STA+ ions on the glass surface is set as u(x, t). The dynamics of x\ and x 2 are written as follows: 



d 2 xi _ dxi dE 



-f3(r ){x 1 (t)-x 2 {t)-r ), (1) 

X — X\ 



d 2 x 2 



dx 2 dE 



dt 2 7 /i0 dt dx 



-PiroKxtW-x^+ro), 



(2) 



where a, ju , and r correspond to the mass of the oil droplet reduced to two lines, a viscous damping coefficient and 
the characteristic size of the oil droplet, respectively. /?(r ) is the elastic constant of the oil droplet with characteristic 
size r (see Appendix A): 

2jw 



ro 



(3) 



E{u) is the surface energy of the glass surface modulated by the STA + ions adsorbed on it: 



E{u) 



eo 



1 + au n 



(4) 



where eo, a and n are positive constants. On the other hand, the dynamics of u can be written as follows: 

du B 2 u 

~dt = dx? + F ( u ' x ' a: iW' a; 2(t);4), 



(5) 



where d is the diffusion constant of the STA+ ions on the glass surface. F(u,x,xi(t),x 2 (t);£ ) corresponds to the 
desorption of STA+ ions from the glass surface to the oil droplet. For this term, we assume the following description: 

F(u,x,x 1 (t),x 2 (t);L ) 

—k\u, i x 2(t) < x < xi(t)) 

k 2 (u a ~u), (£ < x < x 2 (t), xi(t) < x < £- t ) ( 6 ) 
fc 3 (ui -u), (0 < x < £ , I - e < x < I) 
where fci, k 2 , u , I and to are positive constants that correspond to the desorption rate, the adsorption rate, the 
saturated concentration of STA+ ions on the glass surface, the length of the glass substrate, and the length of the 
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FIG. 7: Spatio-temporal plot of the oscillatory motion of a droplet given by numerical calculations using Eqs. <f 1 1 pi - l|18|l . Only 
the parameter Ui was changed: (a) Ui = 0.65 and (b) Ui = 0.658. The other parameters are fi = 0.2, Ro = 0.05, E<j = 5 x 10 -5 , 
A = 1, n = 4, W = 0.1, D = 2.5 x 10~ 5 , L = 0.25, U = 1, K 2 = 1, K s = 1, X w = 0.23, and X 20 = 0.27. The thick lines 
along the x axis correspond to the regions treated with acid. 



glass substrate treated with acid, respectively (see Fig. Eland Appendix B). The initial and boundary conditions are 
set as 



du , „. du . „. 
- M ) = -( M ) = 0, 



(7) 



u(0, x) 



M , (4 < Z < l-l ) 

Ml, (0 < X < i , i-l < X < I) 



(8) 



xi(0) = xxo, x 2 (0) = x 20 (< xxo), 



dxi 



(0) 



dxy 



(0) = 0. 



dt v y di 

To normalize Eqs. - 0, we introduce the following dimensionless variables: 

U=-, X 1 = -, X 2 = -, X=- t , T = k lt , 



(9) 



(10) 
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FIG. 8: Spatio-temporal plot of the slowing of a droplet in the acid-treated region given by numerical calculations using Eqs. 
ill Hi - I|18JL Only the parameter Ui was changed: (a) Ui = 0.68 and (b) U\ = 0.75. The other parameters are the same as 
those in Fig. The thick lines along the x axis correspond to the regions treated with acid. 

Wc then derive the following dimcnsionless system from Eqs. - @): 



dT 2 



dXi dE 
dT~ ~ ~dX 



x=x 1 

B(R )(X 1 (T)-X 2 (T)-R ), 



(11) 



dT 2 



dX 2 dE 



dT dX 



x=x 2 



B(R )(X 2 (T) - Xx(T) + R ), 



(12) 



E(U) = 



En 



1 + AU n ' 



(13) 



dU d 2 U 

Of = D qx2+ F(U,X,X 1 (T),X 2 (T);L ) : 



(14) 
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= < 



F(U,X,X 1 (T),X 2 (T);L ) 

-U, (X 2 (T) < X < X^T)) 

K 2 (U -U), (L Q <X<X 2 (T), 

Xi(T) <X< l-L ) 
K 3 (Ux -U), (0 < X < L , 1 - L < X < 1) 
with the following initial and boundary conditions: 

®L( T 0) = ^-(T 1) = 
dX { ' ; dX [ ' ' 1 



(15) 



(16) 



[7(0, X) 



1, (L < X < 1 - L ) 

C7i, (0 < X < L , 1 -L < X < 1) 



(17) 



X x (0) = X 10 , X 2 (0) = X 20 (< x 10 ), 
« ( 0) = ^(0) = 0, 



dT 



dT 



(18) 



where 



B(Rq 



£ 2 f3(r ) W Txr 2-fw£ 



akf Ro 



, W 



ak\ 



A = aug, 



D = — K 2 = — K = — U = — 

ki£ 2 ' fci ' k\ ' uq 



t - °iV n y - 3:10 y - X20 



VI. NUMERICAL RESULTS 



We performed a numerical simulation using Eqs. (|ll|l - (|15fl under the initial and boundary conditions shown in 
Eqs. OH - (HBJl. 

We only changed the parameter U\, while all of the other parameters are fixed. It is noted that U\ corresponds to 
the saturated concentration at the region treated by acid. Spatio-temporal plots of the droplet given by the numerical 
calculations are shown in Figs. [7| and H3 



12 




WO 



900 
t 



•"j:>: 



lb) 




c 



r 




004 



r. roa. 



;rf< 



MO 



I 1C0 



FIG. 9: Temporal changes in (a) the length X\ — X2 and (b) the velocity of the oil droplet undergoing the oscillatory motion 
shown in Fig. 17Tb) . The high-frequency oscillation is attributed to the eigen- frequency of the mathematical model that did not 
appear in the experimental results. 

When U-y is small enough, i.e., it is difficult for STA + ions to adsorb to the acid-treated region, the oil droplet 
moves back and forth inside the untreated region, as shown in Fig. [7fa). As U\ increases slightly, the oil droplet does 
not turn back quickly but rather rests at the border between the acid-treated region and the untreated region. After 
a brief stop, it begins to move backward, and this intermittent shuttling motion is repeated as shown in Fig. EJb). 
For much larger Ui, the oil droplet does not go back any more, as shown in Fig. |H1 The oil droplet goes into the 
acid-treated region and it continues to move very slowly or stops according to the value of U\ . 

These differences in the characteristics of the motion of an oil droplet are seen only with slight differences in Ui. 
Thus, the numerical results suggest that a slight fluctuation in acid treatment can cause a dramatic change in the 
characteristics of the motion of an oil droplet. In the experiments, the oscillatory motion shown in Fig. |2Jb), the 
stopping of the droplet at the border shown in Fig. I^a), and the slow motion in the acid-treated region shown in 
Fig. EJb) were observed under almost the same conditions. We can guess that the various kinds of motion observed 
in the experiments are due to slight differences in the acid treatment of the glass. 
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We have also noted that intermittent-oscillatory motion (Fig. Efb)) and stopping (Fig. Ufa)) does not appear 
without a diffusion term in Eq. This suggests that the diffusion on a glass surface has a significant effect in 
the real experimental system. There have been several studies concerning the effect of diffusion on the motion of a 
self-running droplet 2, 3^, 3^. Thiele et al. ^(| stated that in their model a running droplet can be stopped by 
increasing the effect of a surfactant diffusion on a substrate. In our study, the diffusion of surfactant on the glass 
substrate is not large enough to stop the motion of a droplet. However, at the boundary of the acid-treated region 
where an oil droplet stopped, the concentration of the surfactant on the acid-treated region is in equilibrium with the 
diffused surfactant from the non-treated region. Thus, the effect of diffusion is optimized by the inhomogeneity in the 
glass surface so that an oil droplet stops its motion. 

Figure El shows (a) the normalized length and (b) the velocity of an oil droplet. When the velocity of the oil droplet 
is large, the length of the droplet increases, and vice versa, which well corresponds to the actual experimental trend 
(Fig.©. 



VII. CONCLUSION 

In the present study, we investigated the self-running motion of an oil droplet on an acid-treated glass substrate. In 
these experiments, the oil droplet undergoes either oscillatory motion within the untreated region or stopping/slowing 
at the border between the acid-treated region and the untreated region. We made a simple "beads-spring" model and 
performed numerical calculations, which suggest that a slight difference in the nature of adsorption in the acid-trcatcd 
region can cause a dramatic change in the characteristics of this motion. The numerical results well reproduced the 
experimental results. In particular, the numerical results suggest that a slight difference in acid treatment can cause 
a dramatic change in the motion of an oil droplet. The results with the model also suggest that diffusion on the glass 
surface plays an important role in the spontaneous motion of a droplet. The inhomogeneity of the glass surface, i.e. 
the boundary of the surface condition, optimizes the effect of diffusion. 

In this study, an oil droplet avoided the area that was treated with acid. It is possible that some other chemicals 
may attract an oil droplet. Above all, the present results showed that a running droplet senses information regarding 
a minute variation in the chemical concentration in space by showing a drastic change in its motion. Thus, this result 
may inspire studies on chemotaxis and chemosensitivity in chemistry as well as biology. 
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FIG. 10: Schematic representation of the framework of the "beads-spring" model. The hemisphere-shaped oil droplet as shown 
in (a) is regarded as a rectangular solid. Since we consider only the x dimension, the conservation law for the volume corresponds 
to the elastic constant of the spring. 
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APPENDIX A: VALIDITY OF THE BEADS-SPRING MODEL 

We discuss here the validity of the beads-spring model for an oil droplet. As shown in Fig. [23, we assume that the 
droplet is shaped as a rectangular solid, whose length, width and height are r, w, and h, respectively. The volume of 
the droplet is conserved: 

V = rwh. (Al) 
The main terms of the energy of the droplet are due to gravity and interfacial tension: 

F =—mgh + ^[rw + 2hw + 2rh] 

oV 2 Q 1 1 

= -^- +1 [rw + 2V (- + -)], (A2) 
Zrw r w 

where 7 is the interfacial tension, p is the density, m is the mass of the droplet, and g is the gravitational acceleration. 
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First, we calculate the fixed point. At the fixed point, dF/dr = 0; therefore, 

'■7[w 5-]=0, (A3) 



dr 2r 2 



which gives the fixed point, tq: 



The elastic constant, /3(fo) is given as 



f P V 2 g 2V\ 1/2 
r ° = I + — ■ ( A4 ) 



d 2 F 



2lW (A5) 



dr 2 

From the above discussion, we can regard the oil droplet as two beads connected with a spring whose elastic constant 
is proportional to the inverse of the characteristic size of the droplet. 

The spring constant is on the same order of magnitude even if the shape of the droplet is regarded as a cap-like or 
ribbon-shaped droplet. 

APPENDIX B: MODEL OF THE ADSORPTION/DESORPTION PROCESS ON A GLASS SURFACE. 

The effect of acids is thought to result in a decrease in the number of adsorption sites X^, r for STA + ions 
Assuming that adsorption-desorption is a langmuir isothermal process, 

STA+ + X- ur ^ STAX sur . (Bl) 

Since mo is the total density of adsorption sites on the glass surface, x is the density of free adsorption sites, and u is 
the concentration of STA + ion in the aqueous phase, the surface concentration it sur of STAX sur under aqueous phase 
changes according to 

^ = ku{x - (1 + ^Hur}, (B2) 

where we assume that x + u SUI = xq, and k and k^ 1 are reaction constants in Eq. (|Bip . The effect of acid is included 
in the decrease in xq, which is equivalent to the decrease in uo to ui in the acid-treated region in Eq. Ip|). 
Similarly, the desorption that occurs under the organic phase is considered as, 

STAX sur + I, -» STAI3 + X- r . (B3) 

We assume that the reaction proceeds only from left to right. Therefore, u SUI changes according to 

= -ku snr v, (B4) 
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where v represents the concentration of Ig in an oil droplet, and k represents the reaction constant. This equation, 
again, is equivalent to Eq. ©. 
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